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Abstract 



' We investigate the phenomenom of emerging regular spectral features from ran- 

r~| . dom interactions. In particular, we address the dominance of L = ground states in 

the context of the vibron model and the interacting boson model. A mean-field anal- 
' ysis links different regions of the parameter space with definite geometric shapes. 

2 ■ The results that are, partly, obtained in closed analytic form, provide a clear and 

transparent interpretation of the distribution of ground state angular momenta as 



observed before in numerical studies. 



1 Introduction 

Recent shell model calculations of even-even nuclei in the sd shell and the pf shell with 
random interactions showed a remarkable statistical preference for ground states with 
L = 0, despite the random nature of the two-body matrix elements, both in sign and in 
relative magnitude A similar preponderance of L = ground states was found in an 
analysis of the Interacting Boson Model (IBM) with random interactions [Q]. In addition, 
in the IBM evidence was found for both vibrational and rotational band structures 0, |^ . 
These are surprising results in the sense that, according to the conventional ideas in 
the field, the occurrence of L = ground states and the existence of vibrational and 
rotational bands are due to very specific forms of the interactions. The studies with 
random interactions show that the class of Hamiltonians that lead to these regular features 

*Sabattical leave 



1 



is much larger than usually thought. These unexpected results have sparked a large 
number of investigations to explain and further explore the properties of random nuclei, 
both for the nuclear shell model [g g g 0, gl, |Tl[ and the IBM g |, |, |12|, 0, . 

In this contribution, we investigate the origin of the dominance of L = ground states 
that emerges from random interactions in two different models: the vibron model and the 
IBM. The vibron model is mathematically simpler than the IBM, but exhibits many of 
the same qualitative features. In this case, most of the results can be obtained in closed 
analytic form. For both models, we suggest to use a mean-field analysis to address the 



problem of the probability distribution of the ground state angular momenta 14 



2 The vibron model 



The vibron model was introduced to describe the rotational and vibrational excitations 

Its building blocks are a dipole boson with = 1~ and a 
The total number of bosons is conserved by the vibron 



of diatomic molecules 
scalar boson with 



m 



0- 



Hamiltonian. We only consider one- and two-body interactions 
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We have scaled Hi by N and H2 by N{N - 1) in 



order to remove the dependence of the matrix elements. The seven parameters of the 
Hamiltonian, e^, e^, Mq, Mi, Cq, C2, fo, altogether denoted by f, are taken as independent 
random numbers on a Gaussian distribution with zero mean. In this way, the interaction 
terms are arbitrary and equally likely to be attractive or repulsive. 

The connection between the vibron Hamiltonian, potential energy surfaces, equilib- 
rium configurations and geometric shapes can be investigated by means of mean-field 
methods lITBIl. For the vibron model, the coherent state can be written as a condensate of 



deformed bosons, which are superpositions of scalar and dipole bosons 
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with < a < 7r/2. The potential energy surface is then given by the expectation value 
of the vibron Hamiltonian of Eqs. (|^) and (Q) in the coherent state 



E(a) 



04 sin"^ a + a2 sin^ a + ao 



(4) 



The coefficients are linear combinations of the parameters of the Hamiltonian 
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For random interaction strengths, we expect the trial wave function of Eq. @ and the 
energy surface of Eq. (^) to provide information on the distribution of shapes that the 
model can acquire. The value of ao that characterizes the equilibrium configuration of the 
potential energy surface only depends on the coefficients 04 and 02. The 0204 plane can 
be divided into different areas according to the three possible equilibrium configurations: 
Si for the s-boson or spherical condensate (ao = 0), 5*2 for the deformed condensate 
(0 < ao < 7r/2), and 5*3 for the p-boson condensate (ao = 7r/2). The distribution of 
shapes for an ensemble of Hamiltonians depends on the joint probability distribution of 
the coefficients 04 and 02 which is given by a bivariate normal distribution 



P(a4,a2) 
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The vectors r and s are defined in Eq. (|5D. The probability that the equilibrium shape 
of an ensemble of Hamiltonians is spherical can be obtained by integrating P(a4, 02) over 
the appropriate range 5*1 (02 > and 04 > —02) 
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16 



0.396 



Similarly, the probability for the occurrence of a deformed shape is obtained by integrating 
P(a4, 02) over S2 (— 204 < 02 < 0) 
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Table 1: Percentages of ground states with L = 0, 1 and A^, obtained in a mean-field analysis 
of the vibron model. 
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Finally, the probability for finding the third solution, a p-boson condensate, is 

P3 = 1 _ _ P2 = 0.388 . (10) 

The rotational spectrum can be obtained from the angular momentum content of the 
condensate in combination with the Thouless-Valatin formula for the moment of inertia 
. The ordering of the rotational energy levels is determined by the sign of the moment 
of inertia 

i?.ot = ^L(L + 1). (11) 

In Table |l| we show the probability distribution of the ground state angular momenta 
as obtained in the mean-field analysis. There is a statistical preference for L = ground 
states. This is largely due to the occurrence of a spherical shape (whose angular momen- 
tum content is just L = 0) for almost 40 % of the cases (see Eq. (|])). The deformed shape 
corresponds to a rotational band with L = 0, 1, . . . , A^, whose ground state has L = for 
positive values of the moment of inertia X > 0, and L = N for X < 0. The third solution, 
the p-boson condensate, has angular momenta L = A^, — 2, . . . , 1 or 0. For X > 0, the 
ground state has L = or L = 1, depending whether the total number of vibrons A^ is 
even or odd, whereas for X < it has the maximum value L = N. The sum of the L = 
and L = 1 percentages is constant. 

Fig. |I| shows that the mean-field results for the percentages of L = and L = 1 
ground states are in excellent agreement with the exact ones. As is clear from the results 
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Figure 1: Percentages of ground states with L = and L = 1 in the vibron model with random 
one- and two-body interactions calculated exactly for 100000 runs (solid lines) and in mean-field 
approximation (dashed lines). 

presented in Table |I], the fluctuations in the percentages of L = and L = 1 ground 
states with are due to the contribution from the p-boson condensate solution. 

3 The interacting boson model 

The interacting boson model (IBM) describes collective excitations in nuclei in terms of 
a system of N interacting bosons [1^. Its building blocks are a quadrupole boson ci^ with 
= 2"*" and a scalar boson with = 0+. The total number of bosons is conserved 
by the IBM Hamiltonian. We consider the one- and two-body Hamiltonian of Eq. ([l|) 
with 




A=0,2,4 
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{d^ xd^Y^^-{dxSY^^ + H.c.] . (12) 



The nine coefficients eg, e^, uq, U2, cq, C2, C4, fo, "^2, are chosen independently from a Gaus- 
sian distribution of random numbers with zero mean |^ . Just as for the vibron model, 
the connection between the IBM, potential energy surfaces, equilibrium configurations 
and geometric shapes, can be studied with mean-field Hartree-Bose techniques by means 
of coherent states [jl^ |1^]. The coherent state can be written as an axially symmetric 
condensate 

I A^, a ) = ^ (cos a + sin a dl^ |0) , (13) 
with —7t/2 < a < 7t/2. The angle a is related to the deformation parameters in the 



intrinsic frame, /3 and 7 [IS, 19|. The potential energy surface is then given by the 



expectation value of the Hamiltonian in the coherent state 

En{ci) = 04 sin^ a + as sin^ a cos a + 02 sin^ a + ao . (14) 

For the IBM, the structure of the energy surface is a bit more complicated than for the 
vibron model due to the presence of the term. This precludes an analytic treatment 
as presented for the vibron model, but qualitatively the results are very similar. In 
practice, for each Hamiltonian the minimum of the energy surface E{a) is determined 
numerically. Again, the equilibrium configurations can be divided into three different 
classes: an s-boson or spherical condensate (ao = 0), a deformed condensate with prolate 
(0 < ao < 7i"/2) or or oblate symmetry (— 7r/2 < ao < 0), and a li-boson condensate 
(ao = ^/2)- Each equilibrium configuration has its own characteristic angular momentum 
content. Even though we do not explicitly project the angular momentum states from the 
coherent state, the angular momentum of the ground state can, to a good approximation, 
be obtained from the rotational structure of the condensate in combination with the 



Thouless-Valatin formula for the corresponding moments of inertia [17|. The results are 
summarized in Table 0. 

(i) The s-boson condensate corresponds to a spherical shape. Whenever such a con- 
densate occurs (in 39.4 % of the cases), the ground state has L = 0. 

(ii) The deformed condensate corresponds to an axially symmetric deformed rotor. 
The ordering of the rotational energy levels L = 0, 2, . . . , 2A^ is determined by the sign of 
the moment of inertia 

^rot = :^L{L + 1) . (15) 

The deformed condensate occurs in 36.8 % of the cases. For T3 > the ground state has 
L = (23.7 %), while for X3 < the ground state has the maximum value of the angular 
momentum L = 2N (13.1 %). 
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Table 2: Percentages of ground states with L = 0, 2 and 2N , obtained in a mean- field analysis 
of the interacting boson model. 
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(iii) The d-boson condensate corresponds to a quadrupole oscillator with quanta. 
Its rotational structure has a more complicated structure than the previous two cases. 
It is characterized by the labels r, and L. The boson seniority r is given by r = 
3nA + A = A^, — 2, . . . , 1 or for odd or even, and the values of the angular momenta 



are L = A, A + 1, . . . , 2A — 2, 2A ||TB[ . In this case, the rotational excitation energies depend 
on two moments of inertia 

Erot = ^r(r + 3) + ^L(L + 1) . (16) 

The (i-boson condensate occurs in 23.8 % of the cases. For X5 > the ground state has 
r = for A^ even or r = 1 for A^ odd (~ 4 %), while for T5 < the ground state has 
the maximum value of the boson seniority r = A^ (~ 19 %). For r = and r = 1 there 
is a single angular momentum state with L = and L = 2, respectively. For the t = N 
multiplet, the angular momentum of the ground state depends on the sign of the moment 
of inertia X3. For X3 > the ground state has L = for N = 3k or L = 2 for A^ 7^ 3fc (~ 9 
%), while for X3 < the ground state has the maximum value of the angular momentum 
L = 2N 10 %). 
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Figure 2: Percentages of ground states with L = and L = 2 in the IBM with random 
one- and two-body interactions calculated exactly for 10000 runs (solid lines) and in mean-field 
approximation (dashed lines). 

In Fig. 1^ we show the percentages of ground states with L = and L = 2 as a 
function of the total number of bosons A^. A comparison of the results of the mean-field 
analysis (dashed lines) and the exact ones (solid lines) shows a good agreement. There 
is a dominance of ground states with L = 0. The oscillations of the L = and L = 2 
percentages with N are due to the contribution of the (i-boson condensate (see Table 
For = 3A; we see an enhancement for L = and a corresponding decrease for L = 2. 
In the mean-field analysis, the sum of the two hardly depends on the number of bosons, 
in agreement with the exact results. 

4 Summary and conclusions 

In this contribution, we have studied the origin of the regular features that were obtained 
before in numerical studies of the IBM with random interactions. In particular, we ad- 
dressed the dominance of L = ground states in the context of the vibron model and the 
IBM. In a mean-field analysis it was found that different regions of the parameter space 
can be associated with definite geometric shapes: a spherical shape, a deformed shape 
and a condensate of dipole (quadrupole) bosons for the vibron model (IBM). For both 
models, we obtained a good description of the probability distribution of ground state 
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angular momenta. 

The studies with random interactions indicate that there is a significantly larger class 
of Hamiltonians that leads to regular, ordered behavior at the low excitation energies than 
was commonly assumed. The fact that these properties are shared by different models, 
seems to exclude an explanation solely in terms of the angular momentum algebra, the 
connectivity of the model space, or the many-body dynamics of the model, as has been 
suggested before. The present analysis points, at least for systems of interacting bosons, 
to a more general phenomenon that does not depend so much on the details of the angular 
momentum coupling, but rather on the occurrence of definite, robust geometric phases 
such as spherical and deformed shapes. For the nuclear shell model the situation is 
less clear. Although a large number of investigations to explain and further explore the 
properties of random nuclei have shed light on various aspects of the original problem, i.e. 
the dominance of L = ground states, in our opinion, no definite answer is yet available, 
and the full implications for nuclear structure physics are still to be clarified. 

In conclusion, the results presented in this article for the vibron model and the IBM 
may help to understand the origin of ordered spectra arising from random interactions, 
as has been observed in numerical studies of nuclear structure models. 
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